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ABSTRACT 

The remanent magnetisation of meteorite material in the solar system indi- 
cates that magnetic fields of several Gauss are present in the protoplanetary 
disk. It is shown that such relatively strong magnetic fields can be gener- 
ated in dusty protoplanetary disks by relative shear motions of the charged 
dust and the neutral gas components. Self-consistent multi-fluid simulations 
show that for typical plasma parameters shear flows with collisional momen- 
tum transfer between the different components of the dusty plasma result in 
a very efficient generation of magnetic fields with strengths of about 0.1-1 
Gauss on spatial scales of two astronomical units in about one year. Based 
on the simulations one may predict that future measurements will reveal the 
strong magnetisation of circumstellar disks around young stellar objects. 
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1 INTRODUCTION 



Meteors are commonly considered to be relics of the protoplanetary nebula in our early solar 
system. Of particular importance in the context of the early evolution of the solar system is 
the remanent magnetisation of meteorite material that indicates the presence of relatively 
strong magnetic fie lds (B = . 



Sugiura et al. 



10 G) in protoplanetary nebula ([Lew &: Sonett 



1978; 



1979). The formation of meteorite chondrules may also demand for magnetic 



fields as the relevant energy source (Sonett 1979; Levy & Araki 1989). Thus, one faces the 
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question of the origin of such strong magnetic fields in protoplanetary disks. Recent numerical 
calculations of collapsing magnetised molecular cloud cores have proven that ambipolar 
diffusion can solve the magnetic flux loss problem of the formation of protostellar objects 



()Ciolek &: Konigll Il993t iDesch fc Mouschoviasl 120011 ) . The magnetic flux can diffuse to the 



region around the pro tostellar object and thus resu lt in a magnetic field of the protoplanetary 



disk of some 0.1 G (IDesch fc Mouschovias 



20011 ) . The fate of this field after the stellar 



object has fully evolved and accretion flow in the disk sets in is, however, not clear. In this 
contribution we consider disk magnetisation in a later stage of the stellar evolution, namely 
after a stellar object with a dipole magnetic field has formed. 

A characteristic feature of partially ionised plasmas is the possibility that magnetic fields 
can be generated without any seed fields. As shown by (Huba & Fedder 1993) a relative 
plasma-neutral gas flow with non-vanishing vorticity results inevitably in the generation of 
magnetic fields. This effect is described by the generalised Ohm's law on the grounds of a 
multi-fluid description. The magnetisation mechanism works significantly more efficiently in 
a dusty plasma (Birk et al. 1996; Birk et al. 2001) where the inertia effects and collisional 
momentum transfer are dominated by the heavy dust grains. 

In this contribution we show by means of multi-fluid simulations that relative shear flows 
lead to a very fast three-dimensional magnetic field self-generation in dusty protoplanetary 
disks. Magnetic field strengths up to one Gauss are generated on the time scale of lyr on 
spatial scales of some astronomical units (AU). 

2 SELF-MAGNETISATION OF PROTOPLANETARY ACCRETION 
DISKS 

We consider the following balance equations that govern the low-frequency macroscopic 
dynamics of a partially ionised dusty plasma (see Birk et al. 1996; Birk et al. 2001) 

-V ■ (p d v d ) (1) 



dp d 
dt 



fh-V-M (2) 

dp n 



dt 

d(p d v d ) 



-V ■ (p n v n ) (3) 
- V ■ (p d v d v d ) - V(p e + pi + p d ) + (pi + Pd)g 



dt 

+ (V x B) x B - v dn pd(v d - v : 



The magnetisation of protoplanetary disks 

- Vinpi (vi - V„) - V en p e (v e - V n ) 

<9(p n v n ) 

— = - V • (p„V n V n ) - Vp n + p n g 

+ VdnPd{v& ~ V n ) + V in pi (v; - V n ) 

+ ^enPe(ve-V n ) 
1 dp e 1 



7e - 1 dt 7e - 1 



V • (p e v e ) - p e V • v e + p n v ne (v e - v n ) 



2 



,PdVde ( k B T e k B T d \ ^piVie ( k B T e k B Ti 



m d V7 e - 1 7«j - 1/ m; \7 e - 1 7* - 1 

PnVne ( k B T e k B T n \ / \2 



2 + pd^de V e - V d 

m n \7e-l 7n-l/ 

PiVi. 

1 dp; 1 



,7e - 1 7n 
+ Pi^e(v e - Vi) 2 



7i - 1 ot 7* - 1 

+ Pifw(vi-Vd) -2 — - 

m d + mi m i + m n \'Yi- 1 7 n - 1/ 



m d + mi 
2- 



m n 


Pi^n(Vi 


-v n ) 2 


+ m n 


Pi^m / 




k B T n 


m, + m n \ 


v7i - 1 


In - 1 









'mj + m d V7i-l 7d-!/ m * W ~ 1 7e - 1 

'Tig , , o 

+ Pi^ieiyi - V e ) 



1 dpd m e . , 2 

r-aT = - rV • (p d v d ) - p d V • v d + Pd^ e (v d - v e ) 

Id ~ 1 at 7d - 1 m d 

+ ; Pd^di(Vd-Vi) H ; Pd^dn(Vd - V„) 

m, + m d m n + m d 



mi + m d \7d - 1 7i - 1 / m n + m d \>y d - 1 7„ - 1 
Pd^de / ^s^d k B T e 



m d \7d - 1 7e 
1 <9p n 1 



7n - 1 dt 7n - 1 



V • (p„v n ) - p n V • v r 



trii 2 m e / 

+ ; Pn^m( V i- V n) H Pn^ne V e -V n 



mi + m„ m 



( k B T n k B T, 

+ ; PnVnd(Vn - V d ) - 2 



- 2 



m n + m d " mi + m n V7„-l 7» - 1 

^ Pn^nd ( k B T n ^ k B T d 
m n + m d V7„-l 7d-l 
' k B T n k B T e 



m n V7n - 1 7e - 1 



4 G.T. Birk, H. Wiechen, A. Kopp and H. Lesch 



rriiC 



e 



V x —\\z d )v id + z d v in )v d - v in v n 

VriiW n d J J 




(10) 



where p a , v a , m a , p a and B are the respective (a = d,i,e,n) mass density, bulk velocity, 
particle mass (the ions are assumed to be protons), thermal pressure (p a = p a kBT a /m a , 
where kn is the Boltzmann constant and T a is the temperature) and magnetic field strength. 
The adiabatic ratios of specific heats are chosen as 7 a = 5/3. The collision frequencies for 
collisions (elastic collisions as well as charge exchange) between species a and (3 are denoted 
by v a p and satisfy p a u a p = ppvpa. By rj the collisional diffusivity (77 = rriiC 2 (u id + u in ) / Anriie 2 ) 
is denoted. Ionisation and recombination are neglected in the present context. The electron 
density is calculated from the quasi-neutrality condition (pi/m, — p e /m e — z d p d /m d = 0) 
where z d is the negative dust charge number. The momentum transfer equations that include 
external gravitational forces with the acceleration g are formulated as the sum of the balance 
equations of all charged components under the assumption of inertialess ions and electrons 
(as compared to the massive dust grains). Most important for the self-generation of magnetic 
fields is the induction equation (see also Ciolek & Mouschovias, 1993) which may be derived 
from the inertialess ion momentum equation. The final term in Eq. (10) is the relevant 
magnetic field self-generation term. The mass, momentum and pressure balance equations 
as well as the induction equation are integrated numerically by the DENISIS code (Schroer 
et al. 1998) under the assumption of ionisation equilibrium. This explicit finite differences 
code is based on a modified Leapfrog algorithm. Extensive tests and applications of the code 
can be found in the papers by Schroer et al. (1998), Birk et al. (2001) and Birk and Wiechen 



Although the present study is motivated by the local information of the magnetised me- 
teorite material in the solar system we assume that our solar system is a standard case for 
normal stellar systems. Thus, in the following we consider a simple protoplanetary accretion 
disk model consisting of electrons, protons, neutrals and charged dust grains using typical 
dynamical parameters for circumstellar disks around T Tauri stars (Montmerle 1991; Konigl 
1994). In particular, we choose the following parameters: n; = 1.01n d = 0.1n n = 10 12 cm~ 3 , 
z d — 1, m d — 100mj, and for the temperatures T n = T { = T d = 500K. The ion- neutral colli- 



sion frequencies can be estimated as (e.g. Krall & Trivelpiece 1986) v in a(Ti)n n ykTi/mi 
where k is the Boltzmann constant and a denotes the cross section. The latter depends on 
the ion temperature in a way that results in a collision frequency that is independent on the 
ion temperature (for more details on the collision theory see, e.g., Mitchner & Kruger 1973). 



(2002). 
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We use a — 5 • 10~ 15 cm 2 as constant approximation of sufficient accuracy in our context. 
The ion-dust collision frequency is v id pz A\f2nn d z1e A \a(\ D / 'a) /3m^ 2 (/cTj) 3 / 2 (Benkadda et 
al. 1996) where \d = \JkT d / %-Kn d z 2 d e 2 denotes the dust Debye screening length. For the 
dust grain size we choose a = 10~ 6 cm. The actual choice of the grain size has no significant 
influence on the induction equation, since a only enters in the logarithm of the formula for 
the ion-dust frequency and does not play an important role for the magnetic field genera- 
tion due to the strong electron depletion (z d — rii/rid ~ in Eq. 10). On the other hand, 
the momentum transfer between the dust and the neutral fluid components depends on the 
grain radius. However, the continuous acceleration (accretion flow) guarantees a finite rela- 
tive dust-neutral shear flow. Without this energy source larger dust grains would result in 
faster saturation and weaker generated magnetic fields. 

Within a region of 10 AU we consider an initially homogenous protoplanetary disk in 
which the dust grains as well as the ions are chosen to be in corotation within a region 
of about 1 AU (Fig.l, upper plot). The corotation of the charged dust grains is caused 
by a stellar dipole field. The neutrals do not couple directly to the magnetic field but via 
momentum transfer with the charged particles. Outside this corotation radius the dust is 
in Keplerian motion (Li 1996). The neutral gas is assumed to be in Keplerian motion with 
the exception of a very inner region of 0.2 AU where the neutral gas is assumed to be in 
corotation to avoid singularities (Fig.l, lower plot). The typical velocities near the transition 
regions between corotation and Keplerian motion are about 1 • 10 6 cm s -1 for the dust and 
1.5-10 7 cm s for the neutral gas, respectively. 

The orders of magnitude of these velocities are in a rough agreement with a central mass 
object of about 2 solar masses, so that we include an external gravitational force correspond- 
ing to such a central mass. The initial configuration is characterised by the immediate onset 
of idealised accretion flow rather than to be in a stationary state. The gravitational field is 
the ultimate energy source for the generation of the magnetic field. 

The simulation box is given by x G [—25, 25], y G [—25, 25] and z G [0, 10] in units of 0.2 
AU. The equidistant grid is resolved by 77 x 77 x 31 grid points. The boundary conditions are 
symmetric at all boundaries, i.e. the values of the scalar quantities and vector components 
at the numerical boundaries correspond to the values at the respective last inner grid points. 

As the upper plot of Fig. 2 indicates, a swirl-like dust flow in the central disk in the 
x-y-plane has built up after 0.7yr. The curl of the relative dust-neutral gas flow (Fig. 2, 
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Figure 1. The initial dust (upper plot) and neutral gas (lower plot) velocity patterns. The velocity field is initially invariant 
in the z-direction, i.e. perpendicular to the disk. The spatial scales are given in units of 0.2 AU. 



lower plot) is sustained by the gravitational pull and results in an ongoing magnetisation of 
the disk. 

The magnetic field after t = 0.7yr and t = 1.4yr is shown in Figs. 3 and 4, respectively. 
The strengths of the different field components are roughly comparable. After t = 0.7yr, the 
amplitude of B z = 0.8G is 5 times the amplitudes of B x and B y and after t = 1.4yr, the 
factor is about 2 (B z = 1.2G). The field generation is caused not only by the poloidal flows 
that develop from the initial configuration but also by vertical relative shear flows arising in 
the dusty accretion disk. 
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Figure 2. The poloidal dust velocity patterns at the middle of the disk (as far as the vertical extension goes) and the z- 
component of the curl of the relative dust-neutral gas flows after ( = 0.7yr. The amplitude is measured in units of the initial 
amplitude of v^ x - 



The magnetic energy density Wb = B 2 /8tt at different heights of the disk are shown in 
Fig. 5. Obviously, the inner part of the disk is completely magnetised in all three dimensions. 

Since the ultimate source of the self-magnetisation is the stellar gravitation the process 
lasts as long as the accretion flow in the circumstellar disk is present. The numerical simu- 
lation, however, fails when the central density gradient becomes too steep (after about 200 
dynamical times). This is a consequence of our simulation approach that does not properly 
deal with the accretion of matter (loss of plasma to the stellar object) and the formation of 
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Figure 3. The magnetic field after t = 0.7yr. The upper plot shows the z-component of the magnetic field and the lower one 
the poloidal magnetic field at z = 5. 

bipolar jets. A global multi-fluid model including all relevant physics is clearly beyond the 
scope of the present state of the art. 



3 DISCUSSION 



We studied the magnetisation of protoplanetary disks that consist of dusty plasmas. Different 
from p revious studies on the magnetisation of the environm ent of collapsing molecular cloud 



cores (jCiolek fc Konigl 



1993 



Desch Mouschoviat 



20011 ) we dwell on the magnetic field 



generation in fully developed accretion disks. The process is driven by accretion, i.e. the 
ultimate source of energy for the magnetic field is the stellar gravitational field, it results in 



The magnetisation of protoplanetary disks 9 





-15 -10 -5 5 10 15 

x 

Figure 4. The magnetic field after t = 1.4yr. The upper plot shows the z-component of the magnetic field and the lower one 
the poloidal magnetic field at z = 5. 

an ordered poloidal magnetic field of some Gauss. It is a very fast and efficient mechanism 
of field generation. Also it is a process of permanent magnetisation as long as some accretion 
flow is present. 

The numerical simulations show that magnetic fields are generated in circumstellar 
dusty disks by macroscopic relative shear motions between the different species. This self- 
magnetisation mechanism that works in addition to the formation of turbulent magnetic 
fields due to the interaction of the stellar field with a turbulent accretion disk (Bardou & 
Heyvaerts 1996). 

We note that in our studies we did not include a stellar magnetic field explicitly. Such a 
field should give rise to a (sheared) wind component close to the star and thus, possibly to 
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Figure 5. The magnetic field energy density (Wb) at different heights (z = 1.4, upper plot; 2 = 5, middle plot; z = 8.6 lower 
plot) of the protoplanetary disk. 
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the further magnetisation of a significant fraction of the magnetosphere. In this contribution 
we concentrated on the generation of magnetic fields in the accretion disks. Further studies 
that include stellar fields are a promising task for the future. 

In our simulations we did not include a gap between the central stellar object and the 
disk for technical reasons. However, since we were interested in the in situ magnetisation of 
the disk material this idealisation does not influence our findings. 

A further idealisation of our study is the initial density homogeneity of the disk. However, 
in the course of the dynamical evolution radial profiles form in all fluid components due to 
the external gravitation. More realistic density profiles perpendicular to the disk would 
demand for the inclusion of self-gravitation. The inclusion of this further energy source that 
also may contribute to the self- magnetisation of the disk is beyond the scope of the present 
investigations but could be of interest for future numerical modelling. Parameter studies, 
however, have shown that the efficiency of the process of magnetisation is not significantly 
influenced by the variation of the initial number densities over a wide parameter range. 

Given the uncertainties of the actual plasma parameters the high efficiency of the dis- 
cussed magnetisation guarantees that protoplanetary disks are expected to generate their 
own magnetic fields of the order of 0.1 - 1 Gauss. In the presence of a small radial mag- 
netic field component the Balbus-Hawley shear instability (Balbus & Hawley 1991) can, in 
principle, result in a faster than linear growth of the poloidal magnetic field component 
provided the appropriate onset criterion is fulfilled. This criterion, however, is not known 
for multifluid systems. On the other hand, localized reconnection may result in a redistri- 
bution of magnetic flux. Investigations on this process are beyond the scope of the present 
contribution and will be considered in further more detailed work. 

If circumstellar disks carry their own significant magnetic fields, dynamic phenomena in 
the coupled star-disk system are rather complex. The disk field should influence characteristic 
phenomena as, e.g., accretion, winds and jets and magnetic activity, significantly. 
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